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Abstract. The mass distributions for giant dipole resonances in ^^S and '^^Sn decaying through 
particle emission and for deep-inelastic collisions between ^^O nuclei have been investigated by 
implementing the Balian-Veneroni variational technique based upon a three-dimensional time- 
dependent Hartree-Fock code with realistic Skyrme interactions. The mass distributions obtained 
have been shown to be significantly larger than the standard TDHF results. 
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INTRODUCTION 

The standard time-dependent mean-field methods used in nuclear physics (namely the 
time-dependent Hartree-Fock (TDHF) approach and its derivatives Q, 0]) have been 
successfully used to determine the expectation values for single-particle observables, 
such as fragment mass, in nuclear reactions and decays but severely underestimate the 
fluctuations in these values [3]. This was first observed by Davies et al. [4] in 1978 who 
performed TDHF calculations to investigate the full-width-half-maximum of the mass 
of the projectile-like fragment for the heavy-ion reactions ^^Kr-i-^^^Pb (£'/a^=494 MeV) 
and ^^Kr-i-^''^Be (£'/a^=600 MeV) and found that they underestimated the experimental 
results by about an order of magnitude. It is now well known that the TDHF method 
underestimates the fluctuation of any single-particle operator, or in fact the expectation 
value of any many-body operator, except the energy [3]. This is a consequence of the 
central assumption of the methods; the approximation of the full many-body state of the 
system as a Slater determinant, evolving in a mean-field, neglecting explicit two-body 
correlations. 

This problem was studied by Balian and Veneroni lH, 0, 0] in the 1980's, who 
produced a general variational theory optimised to the determination of the expectation 
values and fluctuations for arbitrary single-particle operators whilst the state of the 
system is given by a single Slater determinant. They found that, for a system described, 
at the time to, by the one-body density matrix, p (to), the expectation value for the single- 
particle observable Q, at the later time ti, is given by {Q) = Tr [Qp (?i)] , in keeping 
with the usual TDHF approach, whilst the distribution, or standard deviation, AQ, is 
given by 

{^Qbv)\ =rim:^Tr[p{to)-o{to,e)], (1) 



where o e) is a one-body density matrix related to p (t) through the boundary condi- 
tion 

a(?i,e) =exp(ie2)p(fi)exp(-ie2), (2) 

The time evolution of p (?) and a e) is given by the usual TDHF equation. This latter 
result differs from the usual TDHF result 

{AQiDHpf =TT[Qp{ti)Q{l-p{h))], (3) 
h 

in that it depends explicitly on the initial time, with the final time, ti, entering through 
the boundary condition The boundary condition ^ also contains the operator Q such 
that this approach is specifically tuned to the determination of AQ. It has also been shown 
by Das so [8] that, for simple operators such as mass (or, equivalently, charge), the single- 
particle nature of TDHF leads to an unphysical upper limit on the mass distributions that 
can be obtained from ([3]). It can be shown that 

where A is the total number of nucleons in the system. 

Solving ([U) requires that a Hartree-Fock calculation be performed to determine the 
initial state, p (to) followed by a suitable instantaneous excitation of the system and a 
TDHF calculation from to to ti to determine p (?i ) . The transformation Q gives a (?i , e) 
and a second TDHF calculation must then be performed, backwards from ti to tQ, to 
obtain o (to, e). These latter steps must be repeated for a range of values of e to allow 
^Qbv to be determined from the extrapolation of e to 0. 

The large number of calculations required to evaluate ([T]), and their complexity, 
means that only a handful of calculations were performed using this method and those 
calculations used simplified interactions and symmetries (either spherical [|9D, or axial 
[[13, HH]) to make the problems solvable. However, advances in computing power mean 
that this approach can now be implemented using fully three-dimensional TDHF codes 
with full Skyrme interactions EH HE [3]. 



Written in terms of the occupied single-particle states ([U) becomes [llOh 
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[ANbv 

where 



lim^, (5) 



/(e) = E /'dr|(v/„(r,?o,e) |(|)„(r,?o))P, (6) 

m,n<£f 

and A is the number of nucleons in the system. The wavefunctions {f,t,e)) are the 
wavefunctions from the backwards TDHF calculations and are related to the wavefunc- 
tions from the forwards calculations, | ^„ (f , t)), through the boundary condition 

V/(f,ri,e) =exp(ie2)(^)(r,?i). (7) 

In these calculations, as in Bonche and Flocard's earlier work fll'], the wavefunctions 
10/? (^5 ^o)) were obtained by evolving the TDHF equations forwards and then backwards 



without the transformation (|7]) being applied. This approach ensures that all the single- 
particle wavefunctions used in evaluating ([5]) result from the same number of compu- 
tations and has been found to significantly reduce systematic numerical errors [16] al- 
lowing the extrapolation required in Q to be extended to significantly smaller values of 
e. 

We consider the mass distribution in a bounded region of space around a nucleus and 
calculate the mass (number of nucleons) in the nucleus according to 

(A^(i?,))|,= £ [dr\(l>,n{rj)\^0{Rc-\ir-rcM)\): (8) 

m<£f 

where Rc is the cut-off radius used to define the bounded region of space and tcm is the 
centre-of-mass location of the nucleus. 



GDRS IN ^^S AND ^^^SN 

We consider first a giant dipole resonance (GDR) in the deformed nucleus ^^S decaying 
through particle emission. The calculations were all performed in a cubic spatial box of 
size 32 X 32 X 32 fm discretised in steps of 1 fm. The static Hartree-Fock calculation was 
carried out using the SLy6 parametrisation [17] of the Skyrme interaction and produced 
a ground state with prolate deformation = 0. 1 1). A dipole excitation was induced by 
acting on each wavefunction at ? = with a boost 

BD{x,y,z)=Q^^{iFC{A^x+Ayy+A^z)), (9) 

where 




^=\l 1^ / , (10) 

is a spatial cut-off, and where, for protons, F = 1/Z, and for neutrons, F = — 1/(A — Z), 
where A is the nucleus' atomic mass and Z is its charge. The A^, Ay and A^ parameters 
determine the strength of the applied boost and were set to 112.5 fm ^ 

In the present calculation, the system was evolved forward in time until t\ = 250 
fm/c. Dirichlet boundary conditions were used at the edge of the spatial box, where the 
wavefunctions disappear. These lead to spurious reflections for sufficiently large t\ so 
the value chosen for t\ had to be kept sufficiently small, whilst ensuring that (A'^ {Rc))\t 
had stabilised following the prompt de-excitation of the resonance by particle emission. 
The dipole moments, Qx, Qy and Q^, as a function of time, are shown in figure [Tj^a), 



in accordance with [|14(1 



Qi-^^^{{4)-{4)), (11) 

where z = 1, 2, 3 denotes x, y and z and (xf) and (x^) are the expectation values for 
position calculated using the proton and neutron single-particle states respectively. The 
prolate deformation of the ^^S nucleus (with x the long-axis) results in Qy and being 
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FIGURE 1. (a) The dipole moments Qi (i=x, y, z) as a function of time, for a GDR in ^^S. The 
differences between Qx and Qy and are consistent with a prolate deformed ground state where x is 
the long axis, (b) AN^y as a function of e. AA^J^ is given by extraplation to e — > 0. The ;ic-axis has a 
logarithmic scale to emphasise that the values obtained are independent of e across several orders of 
magnitude. The standard TDHF result (calculated at ti and independent of e) is shown for reference. 



identical whilst Qj^ differs. The periodicity of Qx, Qy and Q^ provides an estimate of 
the excitation energy for the oscillations along each of the three primary axes. For Q^, 
a period of ^ 71 fm/c is found, giving an excitation energy ^ 11.5 MeV and, for Qy 
and Qz, a period of ^ 68 fm/c gives an excitation energy Ey ^ 18.3 MeV. 

Following the decay of the GDR in ^^S, it was found that (A^) =26.65 and AA^^^^^ = 
4.08 (using Rc = 8 fm) representing the emission of ^ 5 nucleons. Rc was chosen so 
that the bounded region fully enclosed the nucleus but omitted, as much as possible, 
the radiated components of the wavefunctions. From dH) we obtain {^tdhf) max ~ 
4.46. Following evolution to t\, the transformation (|7]) was applied and the system 
evolved back to Iq for e values down to 10^^. After each calculation the fluctuation, 
AA/|j/(e), was calculated using ([5]) and the results plotted and extrapolated back to 
e = 0. The results are shown in figure [U^b), giving = 5.52 which represents a 

16% increase in AN and exceeds the TDHF upper limit, (AA^|^^^)^^^. Calculations 
were also performed for Rc — 8.5 fm and Rc = 9fm and showed only small changes 
in the results consistent with the bounded region enclosing increasing amounts of the 
wavefunctions tails. 

The calculations were repeated for the doubly magic nucleus ^^^Sn using the same 
spatial box and interaction as in the ^^S calculation. An initial dipole boost dH) with = 
Ay = A, = 600 fm^^ induced the resonance, which was evolved from to = Ofm/c to ti = 
250 fm/c as before. The dipole moments show a periodicity of ~ 88 fm/c corresponding 
to a resonance energy of ~ 14.1 MeV which is close to the experimental value of 
16.1 (7) MeV [[ii. 

The standard THDF calculation gave, at the time (A^) = 121.02 and AA^y^,^^ = 8.46 
representing the emission of 11 nucleons. From dH) we obtain {^tdhf) max ~ l^-^^- 
Again, the Balian-Veneroni transformation was applied for different values of £, the 
extrapolation to £ = yielding AA^J^ = 11.29, representing a 14% increase in AA'^ 
compared with the standard TDHF result. 
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FIGURE 2. Density contour plots showing a deep inelastic collision (Ecm — 128 Me V) between two 
'^O nuclei for an impact parameter b = 2fm(l k, \Ah). 



160+160 COLLISIONS (Ecm = 128 MEV) 

The main application for this method will be heavy-ion collisions as mentioned in 
the introduction. As a test case using light nuclei, collisions between ^^O nuclei at 
Ecm =128 MeV have been studied. 

The dynamic calculations were carried out in a rectangular space of size 50 x 32 x 
50 fm except the head-on calculation where, taking advantage of the symmmetry a 
spatial box of size 50 x 32 x 32 fm was used. The nuclei were initially positioned at 
(±10,0, rkh/l) where b is the impact parameter and the nuclei collide in the x-z plane. 

Calculations were performed for impact parameters b = 0, 1, 2 and 4fm. Contour 
plots showing the densities during a collision for Z? = 2 fm are shown in figure [21 The 
calculations were run for 2000 timesteps (400 fm/c) until the scattered nuclei had clearly 
separated but before they reached the edge of the box. The exception was the calculation 
for b — Aim where the nuclei fused with the compound nucleus decaying through 
particle emission and where the runtime was extended to 2000 fm/c to allow the excited 
compound nucleus time to decay and check that it did not undergo a delayed fission. 

The mass distributions obtained from these calculations using the TDHF and BV 
approaches are given in table [B showing the large increase (at least doubling in each 
case) of /SN in the BV case over TDHF. 



TABLE 1. Mass distributions calculated for ^°0h-'°0 at Ecm = 128 MeV using the standard 
TDHF approach and the Balian-Veneroni variational approach for a range of impact parameters 
(the equivalent angular momentums are also given). 



b(fm) 


1(h) 




ANtdhf 


(ANtdhf) MAX 


ANbv 


Change 





0.00 


14.96 


1.37 


2.82 


2.84 


H-107% 


1 


7.02 


14.90 


1.38 


2.82 


2.87 


H-108% 


2 


14.05 


14.83 


1.44 


2.82 


3.17 


H-120% 


4 


28.10 


25.101. 


2.06 


2.33 


.t 





* After 2000 fm/c. The ""O nuclei fused to form an excited compound nucleus which decays by particle 
emission. The reduced value of {ANt£,hf)max = 4fm is a consequence of the dependence of l|4) on both 
the number of nucleons in the system. A, and on the number of nucleons in the nucleus of interest, {N{Rc))\,. 

^ Not calculated due to the extended runtime required for the compound nucleus to decay and to ensure that the 
compound nucleus did not fission. 



CONCLUSIONS 



The Balian-Veneroni approach has been implemented using a three-dimensional TDHF 
code with the full Skyrme interaction and calculations have been performed for GDR's 
in -^^S and ^■'^Sn and for collisions of ^^O nuclei. The BV approach produces mass 
distributions which are quantitatively larger than those obtained using the usual TDHF 
approach. We are continuing to apply this approach to heavy ion collisions with a 
view towards performing calculations for heavier systems to allow a comparison with 
experimental data. 
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